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Abstract 

We considerably extend our previous analysis of the implementation of an absorptive 
boundary condition, which mimics saturation effects, on the linear CCFM evolution. We 
present detailed results for the evolution of the unintegrated gluon density in the presence 
of saturation and extract the energy dependence of the emerging saturation momentum. 
We show that CCFM and BFKL evolution lead to almost identical predictions after in- 
cluding the effects of gluon saturation and of the running of the coupling. We moreover 
elucidate several important and subtle aspects of the CCFM formalism, such as its rela- 
tion to BFKL, the structure of the angular ordered cascade, and the derivation of more 
inclusive versions of CCFM. We also propose non-leading modifications of the standard 
CCFM evolution which may play an important role for phenomenological studies. 



1. Introduction 

In a previous work, Ref. []]], we have proposed a method for effectively implementing 
saturation and unitarity within a generic linear evolution equation for the unintegrated 
gluon distribution, so like the BFKL [2] and the CCFM 0, [H] equations. The method 
is based on enforcing an absorptive boundary condition at low transverse momenta which 
prevents the gluon phase-space occupation numbers to grow beyond their physical values 
at saturation. Our method is the extension of a strategy originally introduced in relation 
with analytic studies of the BFKL evolution in the presence of saturation 0], whose 
deeper justification lies in the correspondence between high-energy evolution in QCD 
and the reaction-diffusion process in statistical physics [9] . This correspondence is however 
limited to asymptotically high energies and to a fixed coupling [10], whereas our analysis 
in Ref. [l| shows that the absorptive boundary method is in fact more general and also 
very powerful. After reformulating this method in such a way to be suitable for numerical 
simulations, we have demonstrated its efficiency by comparing the numerical solutions to 
the BFKL equation with absorptive boundary condition against those of the proper non- 
linear generalization of BFKL which includes saturation — the Balitsky-Kovchegov (BK) 



equation [11|, We have thus shown that the absorptive boundary method successfully 
reproduces the results of the BK equation for both fixed and running coupling, and for all 
the energies, and not only the asymptotic ones. This success, together with its relative 
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simplicity, makes this method a very compelling tool for phenomenological studies at LHC 
and, in particular, for implementing saturation within Monte Carlo based event generators, 
such as CASCADE pj] and LDCMC [14] . Most importantly, our effective method can 
be also implemented within formalisms whose non-linear generalizations are not known, 
so like the CCFM formalism B a a which lies at the basis of the above mentioned 
generators. It is our main purpose in this paper to provide an extensive numerical study 
of the CCFM evolution supplemented with the absorptive boundary condition, and thus 
demonstrate the role of the saturation effects in that context. Some preliminary results 
in that sense were already presented in Ref. [H , but out present analysis will considerably 
enlarge that previous analysis, in particular by exploring the CCFM formalism in much 
more detail. 

Based on our current theoretical understanding and on extrapolations from the phe- 
nomenology at HERA and RHIC, we expect at LHC a considerably larger phase space 
where saturation effects should be important. The characteristic transverse momentum 
scale for the onset of unitarity corrections is the saturation momentum Qg , and is expected 



to grow quite fast with increasing energy. Next -to-leading order BFKL calculations 15|] 
suggest a power law Q2 ^ gA ^i^-j^ A ~ 0.2 0.3, which appears to be supported by the 
HERA data at small x < 0.01 0, [li]. For forward jet production in proton-proton 
collisions at LHC, one thus expects Qs in the ballpark of 2 to 3 GeV. Even higher values 
could be reached in nucleus-nucleus collisions, or in some rare events, like Mueller-Navelet 
jets [l^. 

A priori, it seems natural to look for saturation effects in the underlying event, that 
is, in the bulk of the particle production at relatively low momenta, where the saturation 
effects (also viewed as multiple scattering [i^, I2H) are clearly important. However, the 
underlying event at LHC will be extremely complex and difficult to study. Besides, for 
such low momenta, it may be difficult to separate saturation physics from the genuinely 
non-perturbative physics in the soft sector of QCD. (A similar ambiguity occurs in the 
interpretation of the small-x data at HERA.) It is therefore useful to recall at this point 
that saturation effects can make themselves felt even at relatively large momenta Q, well 
above Qs, via phenomena like the "geometric scaling" observed at HERA [l^, [ij . [2^. 
Such phenomena, which reflect the change in the unintegrated gluon distribution at high 
k± » Qs due to saturation at low kj_ < Qs Q, 111, S| , are particularly interesting in that 
they represent signatures of saturation in the high-Q^ domain which was traditionally 
believed to fully lie within the realm of the "standard" pQCD formalism — the DGLAP 
evolution 2j] of the parton distributions together with the collinear factorization of the 
hadronic cross-sections. 

Schemes based on the NLO DGLAP evolution have been quite successful at HERA (at 
least for not too low Q^) [25^, but at LHC one will probe much smaller values of Bjorken 
X, and this even for relatively hard Q^ (e.g., in the forward kinematics). The jets to be 
measured at LHC will carry relatively large transverse momenta Q > 10 GeV, but because 
of the high-energy kinematics, their description may require fc^ -factorization 26|] together 
with the BFKL, or CCFM, evolution of the unintegrated gluon distribution. Moreover, 
saturation effects, like geometric scaling, could manifest themselves in hard observables at 
LHC, so like the cross-section for forward jet production 0. Thus LHC will for the ffist 
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time allow us to study saturation physics in the kinematical regime where this physics lies 
in the realm of perturbative QCD. 

The BFKL formalisni, p roperl y generalized to include the non-linear effects responsi- 
ble for gluon saturation [Hf, [H, [23, 0, [3] , is specially tailored to describe the evolution 
of the unintegrated gluon distribution with increasing energy and its approach towards 
saturation. As such, this is well-suited to study the high-energy evolution of inclusive 
cross-sections, and it is able to accommodate important phenomena, like the geometric 
scaling at HERA fl^. [j^ . 23], or the turnover in the DIS structure function F2{x,Q'^) at 
semi-hard 171, However, this is not the right formalism to also describe exclusive 
final states, because it misses the quantum coherence between successive gluon emissions 
in the process of high-energy evolution. Besides, there are additional problems, to be dis- 
cussed at the end of section 12.31 which make the BFKL formalism unsuitable for studying 
exclusive final states. All such problems are corrected in the CCFM formalism [1, Q, 0], 
which has also the advantage to apply within a wider kinematical region, interpolating 
between the high energy evolution (the realm of BFKL) and the evolution with increasing 
virtuality (the realm of DGLAP). A similar formalism derived out of CCFM but using a 
different physical picture for the evolution is the Linked Dipole Chain (LDC) model [30(] 
which covers the same kinematical region as CCFM, and our method of implementing 
saturation can be equally well applied also to this formalism. Infact one of the main equa- 
tions to which we apply our method is equivalent to the master equation in LDC. This 
will be discussed more below. 

Like BFKL, the CCFM formalism is based on the A; _l -factorization. This makes it 
possible to take into account some of the NLO corrections in the collinear approach by 



simply treating the kinematics of the scattering more accurately 3l|, l32l . l33l |. Still like 



BFKL, the CCFM evolution resums all powers of ag which are accompanied by large energy 
logarithms Ins, with s the center of mass energy. Li fact, the CCFM and BFKL evolutions 
yield identical predictions for the dominant behaviour in the formal high-energy limit s — > 
00. But this formal limit is conceptually unrealistic and phenomenologically irrelevant, 
since it violates unitarity. What is relevant, is the approach towards the unitarity limit and 
gluon saturation with increasing energy, which is a priori different in the two formalisms. 
In the recent years, this approach has been extensively studied within the BFKL evolution. 



by using its non-linear generalizations: the BK ll|, [l^ and the JIMWLK equations 



271 . l28l |. However, no such a study was performed within the context of the CCFM 



evolution prior to our recent work [l,]. As already stated, such a study is the main objective 
of the present work. 

The present study should be viewed as a step towards a systematic inclusion of the 
effects of saturation in the description of exclusive final states. Most studies have so far 
concentrated on more inclusive observables, for which a good description can generally be 
obtained (for sufficiently high Q^) also with linear evolution equations alone. A number 
of papers emphasized the need and importance of studying saturation in exclusive final 



states [3J, |35|, |36|, |37|, but so far the explicit studies were mostly confined to inclusive 



observables. The possibility of looking at more exclusive observables will undoubtedly 
make it easier to distinguish the predictions of linear and non-linear evolutions. 

But before we consider the effects of saturation, we shall dedicate a large part of this 



3 



paper to a detailed presentation of the respective linear evolution, in order to clarify several 
non-trivial aspects of the CCFM formalism which are important for our purposes. In doing 
so, we will try to carefully motivate the steps leading to the final evolution equations 
(simplified versions of the CCFM formalism) on which we shall apply our saturation 
boundary. Among the different aspects of CCFM to be discussed here, there are parts 
which have been already presented in previous papers that we shall refer to, but there are 
also parts which to our knowledge have never been presented before. In this paper, we shall 
try to give a unified presentation of all the relevant aspects, using a intuitive geometrical 
representation for the phase-space of the CCFM evolution. To facilitate the reading of 
the paper, we have moved some of the most technical developments to appendices, and 
kept only the important results in the main text. Here is a summary of the topics to be 
covered in what follows and also of our main conclusions: 

Sect. 2 will introduce the basics of the CCFM formalism, that the rest of the work will 
rely on. In particular, in Sect. [2^2] we shall clarify the relation between the phase-space of 
the CCFM evolution and that of the BFKL evolution, thus recovering previous results in 
Ref. 381], but from a different perspective. In Sect. l2.3l we shall present the standard version 
of the CCFM evolution as an integral equation, and on this occasion we shall explain the 
approximations which are involved in this rewriting and which are often kept implicit in 
the literature. In Sect. 12.41 and Appendix A we discuss some subtle aspects, and correct 
some mistakes in the literature, concerning the virtual ('non-Sudakov') form factors which 
express the probability for not emitting gluon in the course of the evolution. The careful 
derivation of these form factors, as outlined in Sects [2?3l and [2741 and in Appendix A, will 
also allow us to better understand their physical origin and thus propose some improved 
expressions for them, which treat more accurately the kinematics (by including effects 
of recoils in the energy). The new form factors, to be presented in appendix B, differ 
from the standard ones by terms which are formally of higher order, but which may be 
numerically importan10 and thus interesting for the phenomenology. Another aspect of 
the form factor is its appropriate form in the formal high energy limit, which is relevant 
again for the comparison to BFKL and this will be discussed in appendix C. 

Sections [3] and m are the key sections in this paper. In Sect [3] we successively simplify 
the CCFM formalism and reduce it to a set of simpler, integral and differential, equations, 
which are more suitable for numerical simulations. The most general equation, the integral 
equation ()3.ip . explicitely preserves all the hard and soft gluon emissions from the t- 
channel propagators. This is the equation at the basis of the CASCADE Monte Carlo event 
generator [l3j], and it can be generalized to include saturation effects, as we shall explain in 
Sect. 14.11 But for the present purposes, it is preferable to work with simpler versions of the 
CCFM evolution, namely Eqs. p. lip and (j3.12p . which are more 'inclusive' — in the sense 
that some of the virtual form factors are used to cancel the 'soft' gluon emissions. Such 
cancellations are not exact, but rather require some additional kinematical approximations, 
which are however in the spirit of the CCFM formalism. These approximations are also 



^We shall not include these new form factors in our present numerical analysis since their structure is 
such that they can be efficiently implemented only within Monte Carlo simulations. 
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similar to those underlying the 'Linked Dipole Chain' (LDC) model |3QI . And indeed, 
our Eq. (j3.1ip appears to be equivalent to the master equation in Ref. [30|], although our 
respective approaches are quite different. Intriguingly, it turns out that if one performs 
the same type of approximations starting with the BFKL formalism, one is again led to 
the same two equations (I3.11|) and (I3.12|) (see the discussion in Sect. 13.41) . In our opinion, 
this points out towards a deep similarity between the BFKL and CCFM formalisms, which 
when applied to inclusive observables differ only in the accuracy to which they treat the 
kinematics. 

In Sect. [4] we introduce the absorptive boundary method which effectively implements 
saturation within a linear evolution equation, so like BFKL or CCFM. First, in Sect. 14. 1^ 
we describe this method on the example of the BFKL equation, where the comparison 
with the non-linear BK equation will allow us to demonstrate the success of the method. 
Then, in Sect. 14.21 we analytically study the high-energy behaviour of the approximate 
CCFM equations ()3.1ip and ()3.12p and thus determine the energy-dependence of the 
associated saturation momenta, in the case of a fixed coupling. Our analysis shows that 
the CCFM evolution is somewhat faster than the BFKL one: the respective saturation 
exponent is slightly larger. On the other hand, the characteristic functions which determine 
the momentum-dependent anomalous dimension, are very similar to the BFKL one, for 
both Eq. (j3.1ip and Eq. (j3.12p . In Appendix D, we relax some of the approximations 
used in deriving Eq. ()3.1ip and thus obtain a more accurate equation, whose high-energy 
predictions are even closer to those of the BFKL equation. 

Finally, Section [5] presents a systematic numerical analysis of the various evolution 
equations — the BFKL equation and the simplified versions, Eqs. p. lip and (j3.12p . of 
the CCFM evolution — with the purpose to illustrate the role of the saturation boundary 
and also of the running coupling effects. We first demonstrate that the respective linear 
equations (no saturation boundary) lead to rather different evolutions, which are moreover 
infrared unstable in the case of a running coupling. Then we show that the infrared 
instability is cured after including the saturation boundary (the saturation momentum 
effectively acts as a hard infrared cutoff which increases with the energy) and moreover 
the respective predictions of the various equations remain very close to each other, up to 
astronomically high energies. Hence, in so far as the unintegrated gluon distribution is 
concerned, the BFKL and CCFM evolutions properly corrected for saturation and endowed 
with a running coupling are rather similar to each other. 

2. Review of CCFM 

Our aim in this section is to review the CCFM approach, mainly the work in This 
will prepare us in understanding our subsequent strategy for simplifying this formalism 
and most efficiently complete it with a saturation boundary. The original formulation 
in [3] is very careful and complete, but also quite technical and not always transparent. 
We shall therefore try to mostly give a geometrical representation of the equations to be 
presented here. In this process, we will derive to some interesting results that we were not 
aware of, and also clarify some points which are often confusing in the literature, such as 
the precise form of the "non-Sudakov" form factor A^^. 
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Figure 1: Kinematics of the gluon radiation. Tie last t-channel propagator is denoted by Q„, wliile the 
real s-channel gluons are represented by the horizontal lines. 

2.1. Kinematics and Basics 

We use Fig. [1] to define the kinematics and schematically introduce the physical pic- 
ture. This figure represents a gluon ladder as produced by the CCFM evolution; we denote 
by qi the transverse momenta of the real, s-channel, gluons, and with Qi the transverse 
momenta of the space-like, t-channel, propagators, which are not explicitly shown in the 
figure. The incoming virtual gluon has zero transverse momentum Qo = (that is, this 
gluon is taken to be collinear with the parent hadron, not shown here), and hence one has 
Qi = — X]fc=i Qk- We will use yi and Xi to denote the energy fractions of the s-channel 
and t-channel gluons, respectively, measured with respect to the energy E of the incom- 
ing proton. In Fig. [H the s-channel (or 'real') gluons are enumerated according to their 
energy: 

X = Xn <. Un <. Vn-l <C • • • < 2/2 < 2/1 ~ 1 , (2.1) 

but this ordering is not necessarily the same as that of the gluon emissions along the ladder 
(i.e., it is not assumed that the gluon with energy fraction yi is emitted out right after 
that with fraction yi-i, etc.). Rather, as we shall shortly see, the real gluon emissions in 
the CCFM ladder are ordered according to their angles = qf/{yfE'^). This ordering 
issue is potentially confusing, since e.g. the relation Qi = — qu would be strictly true 
if the labels i attached to gluons in Fig. [T] were also indicating their order of emission, 
i.e., if the real gluons emissions were ordered according to their energy (which they are 
not). The resolution of this puzzle, to be explained in detail later on, is that the actual 



emissions which do not obey energy ordering are also soft in the sense of carrying httle 
transverse momenta -conserving'), and hence they do not affect the momenta Qi of 
the i~channel gluons: the latter are fully determined by the 'hard' emissions which are 
simultaneously ordered in angle and in energy. (See Sect. 12.31 for details.) 
The (integrated) CCFM gluon structure function can be written as 

n i=\ \ ^* perm / 

^—b{x - x^)Sl,^{y^~i)Sl,{Vl ...n) (2.2) 

where yn,+i = 0, and S^ik and Sm are the virtual corrections associated with the eikonal 
and the non-eikonal vertices respectively. The theta function 0(^/, is a consequence 

of the quantum coherence between successive emission which implies that the emission 
angle must increase when moving upwards along the ladder (i.e., towards the hard scat- 
tering). Notice that one can have any angular ordering for the given energy ordering. 
(We shall later relabel the real gluons according to the angular ordering.) The angle ^ 
in the argument of A is the maximum angle allowed by coherence and is determined by 
the kinematics of the hard scattering; roughly, ^ ~ Q"^ /x^E'^ with the virtuality of the 
incoming photon. In this case, the structure function (j2.2j) gives the gluon distribution, 
A{x,i)=xg{x,Q^) [39]. 

The virtual form factors in (j2.2p are given by 



Seik{yi,0 = exp --aW — / -r\ (2-3) 




Sne{l2...n) = Ue^pilas r — f ^ ) ^ TT 5ne(A;) (2.4) 



where ^{Qk) = Q\/{y'^E'^) and the integral over the transverse momentum cf' has been 
written in terms of the associated angular variable = (f' /y'^E'^ . Notice that the exponent 
in the non-eikonal form factor is positive. The apparent divergence in Seik aud in the real 
emission density d^/^ is regulated by a collinear momentum cut go (which also regulates 
the dy/y divergence in Seik)- 

One should be aware that (j2.2p does not correspond to an exclusive final state. In 
p.2p one has already inclusively summed over all subsequent emissions from the outgoing 
gluons Qi. Each such gluon can further radiate within a cone of half opening angle 
These final emissions are such that the real emission probability is exactly compensated 
by virtual corrections of the type S^ik , and they are therefore not visible in the expression 
(12. 2p for the (inclusive) gluon distribution. For the study of exclusive final states, however, 
one needs to include all the emissions. 

2.2. The phase space and relation to BFKL 

In this subsection we shall describe in detail the phase-space for real and virtual gluon 
emissions within the CCFM ladder and then argue that, up to subleading effects, this is 
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Figure 2: Geometrical representation of two emissions (fat dots) in CCFM shown in the (ln(l/t/), In q^)- 
phase space. The total phase space is bounded by the diagonal line which indicates the maximal angle 
allowed by coherence. The shaded regions show the phase regions over which Sne{k) (left) and Seik (right) 
are integrated over. 



essentially the same as the phase-space for the BFKL evolution. This will allow us to 
conclude that the CCFM and BFKL evolutions become identical with each other in the 



high energy limit. (See also Ref. 38| for similar considerations.) 

Let us start with the virtual emissions and rewrite the eikonal form factor in (12.31) as 



Seik{yi,i) 




(2.5) 



k=l 



with yn+i = 0. Let us now consider the two emissions in Fig. El In this figure the horizontal 
axi^ is ln(l/y) while the vertical axis is ln(g/gg) (from now on we shall for simplicity 
omit using qq). A black dot in such figures denotes a real gluon with the respective 
values for the energy (y) and the transverse momentum {q'^). Since In^ = ln(g^/y^i?^) = 
ln(g^/£'^) + 21n(l/y), the emission angle will be constant along diagonal lines in the figure. 
The diagonal line shown in figure [2] denotes the maximum angle, determined by ^. The 
diagonal lines parallel to this line and which pass through the gluons will indicate the 
angle of the gluons (see e.g. Fig. [5]). 

The shaded regions in the two figures indicate the phase space over which the non- 
eikonal and eikonal form factors are integrated over. The rightmost figure, representing 
Seik{k), is easier to understand. Here the y integral is bounded by yk+i and y^, that is we 
integrate over the region between the two vertical lines, and the ^ integral goes up to ^, as 
can be seen in the figure. To understand the leftmost figure, note that in Sne{k) defined 
in (12. 4p . the y-integral is the same as in Seik{k) while the .^-integral (or, equivalently, that 
over = y^i?^^^) is integrated from Qf. up. The horizontal line bounding the shaded 
region in that figure corresponds to InQ^. Since the exponent is negative in Seik{k) and is 



These type of diagrams have been widely used for example in |30t|, but notice that there the horizontal 
is taken as In(^) instead of ln{l/y). Our choice is the same as the one in [3^ . 
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Figure 3: The shaded region is the region left over where there is a net contribution from the virtual form 
factors. It is bounded from above by Qk- The real emissions are represented by fat dots. 



positive in Sne{k), we get a complete cancelation in the region of overlap. It is then easy 
to see from Fig. [2] that the region which is left is the one shown in Fig. [3l and that the 
net exponent is negative. 

Now, remember that in the BFKL evolution, the virtual corrections are contained in 
the "non-Sudakov" (or "non-eikonal" ) factor A^f^^{k) defined by 

and we see that this corresponds to the shaded area in Fig. [3l Thus we find 

Sl{k)-Si,{k) = A^r^{k). (2.7) 

Consider now a complete set of emissions in the initial chain as shown in Fig. HI In 
the figure the gluons are obviously enumerated according to their energy. The explanation 
of the phase space in figure is the following. The Qi are determined by the relation 
Qi = — Ik ■ For example, since Qq = 0, we have Qi = qi (in the following all the 
momenta denote the norm of the transverse components). Since g2 S> qi we have Q2 — 92; 
where the small recoil is neglected in the figure. Then the subsequent real gluons have 
small momenta so that Qk+i ~ Qk up to gluon 6 which has large momentum, and therefore 
Qq ~ qQ and so on. In the end we get an integral over the total shaded region in Fig. HI 
and again this is exactly the same that we would have in BFKL. 

Of course in CCFM the total phase space is determined by ^ so in Fig. U] we have 
assumed the constraint ^ < ^ not to cut the shaded regions. Hence, a difference between 
the CCFM and BFKL ladder appears towards the end of the chain. This difference however 
is not enhanced by ln(l/a:), and hence it is subleading from the viewpoint of the BFKL 
resummation. Moreover, the CCFM real-gluon emissions are ordered according to their 
angle, and not to the energy. However, in the ensuing gluon distribution (j2.2p there is 
a sum over all the possible angular orderings for a given energy ordering. Therefore the 
phase space for real emissions is also the same, up to the subleading difference mentioned 
above. 
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In general, however, the energy weighting of the real and virtual emissions is different 
for the BFKL and CCFM evolutions. For the latter, this is encoded in the splitting 
functions in Eqs. ()2.2p and ()2.3p - ()2.4p . which show that the (real and virtual) gluon 
emissions are distributed logarithmically in y — the rapidity of the s-channel gluons. The 
BFKL evolution, on the other hand, retains those diagrams which resum all orders in 
aslji{l/x) : these are gluon ladders in which the t-channel gluons are strongly ordered 
in longitudinal momentum {x = Xn <^ Xn~i ^ • • • ^ xi <^ xq) and distributed with 
the logarithmic weight dxi/xi. Clearly, strong ordering in the t-channel implies a similar 
ordering in the s-channel — from Xj <C a^i-i, it follows that yi ~ <C yi-i ^ Xi-2 
— , so that the positions (in energy) of the t-channel propagators uniquely determine the 
positions of the real gluons. This is true in the strict high energy limit, where s — > oo 
at fixed Q^, and Os is very small. Beyond this, however, it is important that the CCFM 
configurations are generated according to the rule of quantum coherence, and hence they 
represent realistic final states (at least, up to further emissions from the real CCFM gluon, 
as explained as the end of Sect. 12. ip . Thus, although the two types of evolution provide 
identical results for the (inclusive) gluon distribution in the formal high-energy limit the 
CCFM evolution is more appropriate for describing actual final states (see the discussion 
at the end of section [2?3l) . Moreover, this formal high-energy limit becomes meaningless 
in the presence of saturation, as we shall later explain, and when this limit is properly 
taken (see in Sect. S]), differences are expected to appear already in inclusive quantities, 
so like the gluon distribution. 

2. 3. The structure of the angular ordered cascade 

Returning to Eq. (j2.2p . this can be further simplified to obtain a more familiar expres- 
sion for the gluon distribution. To that aim, we shall divide the initial state radiation into 
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Figure 5: Dividing the CCFM radiation into two subsets: hard emissions (enumerated emissions) and soft 
emissions (indexed by small letters) as explained in the text. All real emissions are represented by fat dots. 

two classes, "soft" and "hard" (or "fast") gluons [1]. This is done as fohows. Consider 
the set of initial gluons shown in Fig. [D The "hard" subset of gluons are those which are 
not in angle followed by a gluon with more energy, i.e. with higher y. All other gluons 
are defined as being "soft" gluons. In Fig. [5l the gluons marked by 1, 2, 3, 4, 5 and 6 
are hard gluons, since as compared to them, there are no other gluons with larger angle 
and higher energy. (Recall that a larger angle would mean a gluon above the respective 
diagonal line, while a higher energy means a gluon to the left of the vertical line through 
the gluon.) The gluons marked by a, b, c, d, e, /, g, h, i and j are soft gluons. For example, 
a is in angle followed by 1 and 1 has larger energy than a since it is located to the left of 
it. Gluon d is followed by e which has larger energy, which itself is followed by 4 in angle 
which has also larger energy than e. The hard gluons are thus ordered both in angle and 
in energy. 

The soft gluons can furthermore be divided into clusters. Define the cluster Ck as 
consisting of those soft gluons which have their angle between (,k-i and and which 
have energies less than y^. Thus in Fig. [H a belongs to cluster Ci since < < '^i 
(^0 = 0) and Va < Hi, b and c belong to C2 since ^1 < (,b,c < ?2 and so on. In this example 
there are no gluons in C3. The phase spaces for the clusters Ci are shown in Fig. [HI We 
now see that the "soft" gluons which belong to the cluster Ck are indeed soft in the sense 
of having lower energies than the "hard" gluon k which defines the cluster. 

The advantage of this separation is that it allows us to rewrite the distribution ()2.2p in 
a simpler way. To that aim, notice that the soft gluons in have all transverse momenta 
smaller than g^. This is obvious from the figures and more formally follows from the fact 



11 



log q"^ 




Figure 6: The representation of the hard emissions, enumerated fat dots, together with the clusters d. 
For each hard emission i, the cluster d contains all the soft emissions associated with that hard emission. 

thalU ^s,k < Cfc and ys,k < Vk implies ql^^ < {ylk/Vkhk < Qk- Then one can write 

SL = n ^neik) « n Slik) (2.8) 



where £/ denotes the set of all emissions while denotes the subset of hard emissions, 
and the phase-space for a hard gluon in Jif is the whole respective cluster Ck, as shown 
in Fig. [6j The gluon distribution (|2.2p can then be written as 



E/n 

n=l-^ k=l 



as — ^Sl^{yk+i,yk,Qk)diCk+i - ik)6{Vk - Vk+i)] —5{x - x„) 
Vk ik J Xn 



n+l / oo r ^ A AC \ 

X SUy.rOWiY.^"^ W—^^e{i,.^,-^A. (2.9 

k=i \m=o ■'^h i=i J 

where ^n+i = aiid we now index the gluons by using the angular ordering. The meaning 
of this equation is simple: it says that we can construct each chain by adding an arbi- 
trary number of soft gluons between each pair of hard emissions. We can now further 
simplify this expression, by using cancelations between real and virtual contributions to 
the emission of the soft gluons. 

To that purpose, we refer to Fig. [7]where we have identified the hard and soft emissions 
from Fig. [SJ and the numerated emissions are the hard ones. For these we also indicate 
the angles. The soft gluons are marked by small letters. We now define the new "non- 



^We denote by ^s,fe and ys^k the angle and energy fraction for a soft gluon belonging to cluster Ck. 
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Sudakov" form factor A„s by 

Ans = Y[Ans{k) =l[exp{-asAk). (2.10) 

k k 

Similarly we define the "Sudakov" factors As{k) by 

A, = []A,(A:) = J]exp(-Q,Cfc). (2.11) 

k k 

Thus we have 

Si,{yi,0 - U^neik) = llAUk) . A,{k). (2.12) 

k k 

Now, the summation over the real soft emissions in each cluster in (j2.9p exponentiates, 

m=0 •^'^fe j=l 

(2.13) 

By the definition of the Sudakov in (|2.11|) we thus have 

A,{k) . A'j^f\k) = 1, (2.14) 

i.e. the real soft emissions are exactly compensated by the Sudakov form factors. After 
thus inclusively summing over all soft emissions, the structure function can be finally 
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written as 



oo I. n / 

^(^'0=E/ n 

n=\-' k=\ ^ 



Vk ik 



\Ls{yk+i,yk,Qk)0{£,k+i - CkWivk - Vk+i) — 5{x - Xn) 



(2.15) 



This expression involves the hard gluons alone. 

Before leaving this section, two comments are in order: 

(i) In Fig. [7] the horizontal lines between the real emissions represent the momenta 
of the t-channel gluons. As obvious from the figure, these momenta are determined solely 
by the hard subset of emissions (c./. Eq. (|2.8p ) which are ordered in both energy and 
angle. This a posteriori explains why the condition Qi = — qu is approximately true 
irrespective whether the labels i, k refer to energy ordering, or to the angular ordering 
(which is the actual order of the gluon emissions). This argument shows that there is an 
implicit approximation in the CCFM formalism — the fact that soft gluon emissions are 
assumed not to change the virtual transverse momenta. Hence, without further loss of 
accuracy, we will later use similar approximations to simplify the expression of the gluon 
distribution even further (see Sect.Ej). 

The second important point is the relation to BFKL mentioned earlier and which 
deserves some clarifications. To compare to BFKL, it is convenient to define Zk by Xk = 
Zj^Xk-i- This implies = X}^_i[l — z^)), and therefore 



Thus in the CCFM ladder one can distinguish two vertices contributing to the split- 
ting Qk-i Qk + Qk- The one corresponding to the small-z^ pole 1/zk is dubbed the 
"non-eikonal" vertex as it comes from the piece of the three gluon vertex in which the 
polarization of the parent gluon is inherited by the real gluon q^. The opposite 1/(1 — z^) 
pole is dubbed the "eikonal" vertex, and in this case the polarization, together with most 
of the energy, is inherited by the t-channel propagator Q^. The hard emissions previously 
identified are associated with the l/z/. pole in Eq. (|2.16|) , while the soft emissions with the 
1/(1 — Zk) pole In contrast, in a BFKL ladder, only the l/z^ pole would be present, 
and the energy ordering coincides with the actual sequence of emissions along the cascade. 
In that case the typical gluons are such that <C 1 and hence ~ Xk-i, as anticipated at 
the end of Sect. 12.21 Thus in the corresponding phase-space integrals, like Eq. (12. 6p . one 
can replace the measure dyk/vk by dz^/zk- We are now prepared for our second comment: 
(ii) Despite the formal equivalence between the CCFM and BFKL evolutions (in the 
high energy limit), the latter cannot be used to generate the final state, not only because 
it does not obey the condition of angular ordering (as required by quantum coherence), 
but also because the Regge kinematics z^ <^ 1 cannot be ensured in practice when trying 
to generate a BFKL ladder. The reason is as follows: the BFKL emission probabilities for 
real and, respectively, virtual gluons are separately infrared divergent (see e.g. Eq. (12. 6p ) 
and thus require an infrared regulator qQ. Although the dependence upon go formally 
cancels in the complete result, the introduction of this soft momentum cutoff will falsify 




(2.16) 
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the condition that ^ 1 in the intermediate steps. Indeed, when a new value for is 
randomly generated with probability law Ail^^^^\k), the typical value value is such that 

In- ^ I , (2.17) 

Zk as ln{Ql/q^) 

which in principle should be of ©(l/a^) for the Regge kinematics to apply, but in reality 
becomes of 0(1) (meaning Zk ~ C(l) as well) whenever qq is taken to be small enough. 

Within the CCFM evolution, this problem is avoided due to the presence of both 
types of poles, 1/zk and 1/(1 — Zfc), and to the angular ordering. In that context, the 
dependence on qq is present in the Sudakov form factor, as the areas in Fig. [7|are cut 
from below by qq. In the qo ^ limit, successive emissions will become very close to each 
other in angle. Indeed, an emission typically occurs when the corresponding region 
(that we now define separately for each emission, either hard, or soft; see Fig. [8]) has an 
area of 0(1), 

/ «2 [J, \ r~c. 

- 1, (2.18) 

which implies ~^ (,k-i when qq — > 0. In that case one can identify two limiting cases: 
(i) z <^ 1 and q more or less similar to the momentum of the previous emission, or (ii) 
2 ~ 1 and q either of the order of, or much smaller than, the q of the previous emission. 
These two possibilities are precisely the type of emissions already present in CCFM and 
therefore the structure of the cascade is not altered by unlike what happens in BFKL. 

Of course, in the presence of saturation the dependence on any soft momentum disap- 
pears naturally, as the dynamically generated saturation momentum, which grows rapidly 
in the course of the evolution, provides a natural cutoff (see the discussion in Sect. 14. ip . 

2.4- The virtual form factors 

In this subsection we shall display some more explicit formula for the non-Sudakov 
and Sudakov form factors Ans and A^, whose detailed derivation is presented in appendix 
A. Although such formula were already presented in the original work 0], it turns out 
that they are often written in a wrong way in the literature (see the discussion in appendix 
A). To avoid such errors, and also in order to be able to generalize these form factors by 
including non-leading effects — a task that we address in appendix B — , it is essential 
to have a proper understanding of the derivation of the corresponding formute. This is 
briefly discussed here and then in more detail in appendix A. 

Let us first recall that for the hard emissions, = (1 — Zk)xk-i ~ Xk-i, while for the 
soft emissions x/. = z^Xk-i ~ a^fe-i- Then one can as a first approximation set = for 
the hard emissions (except in the l/z^ pole), and z^ = 1 for the soft emissions (except in 
the 1/(1 — Zk) pole). Therefore the region A^. has the transverse momentum bounded from 
below by and from above by Qk, while the respective y integral is bounded between 
yk and Vk+i- Thus we approximately have 

A„,(fc)=exp -a, / ^ -f). (2.19) 



15 




Figure 8: Hard and soft emissions, with the Sudakov (regions d) defined for each individual emission. To 
leading order there is no non-Sudakov (regions Ai ) associated with the soft emissions. Real emissions are 
again represented by fat dots, while the k± and energy fraction of the virtual propagators are denoted by 
crosses. In this example, emission number 2 is a soft emission. 



Defining y = zxk-i, and switching from ^ to g by using ^ = q /{y E ) one gets 

^ns{h) = exp {-as C — ^1 . (2.20) 

The reason we did not set 1 — Zfc = 1 in the lower hmit of the q integral is because the 
CCFM equation is usually written in terms of the so-called rescaled momenta defined by 
Pk = qk/{^ — Zk) so that the factor (1 — Zk) is absorbed into the definition of pk- Then 



A„,(fc) = exp {-as — ^ (2.21) 




which is the form used in [J]. Equation 12.211 for A„s is, however, usually written in a 
different way in the literature (see Eq. (jA.ip in the appendix). In appendix A, we give 
a more careful derivation of the non-Sudakov form factor, and we demonstrate that the 
correct form is indeed given by 12.211 

To write down the Sudakov one first needs to define it for each, hard and soft, individ- 
ual emission. In Fig. [8]we show an explicit chain of hard and soft emissions, where the real 
emissions are indexed according to their angular ordering, and where we also explicitely 
show the virtual t-channel propagators by crosses. The individual Sudakov form factors 
are then defined as the integrals over the regions Ck in the figure (note that these are not 
the same regions Ck as before, now we define such a region for each emission, not just for 
the hard ones). 
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We see that the region Ctis to the left bounded by the energy of Qk-i, while there 
is no lower limit for the energ}0. In momenta it is bounded between the angles of the real 
gluons Qk-i and qk- Therefore we may write the Sudakov as 



where e represents the soft cutoff. Since we must have g > go we get ^ = q^/{'ip'E'^) > 
qQ/{y'^E'^) and y > qo/{VCE) = e. Defining p = y/^Xk-iE, we have e = Xk~iqo/p and 



where e' = e/xk-i = qo/p- Now if one wishes one can let y = 1 — z and then the usual 
form for the Sudakov form factor is obtained. 

3. More inclusive versions of CCFM 

In section 12.31 we have shown that one can use the Sudakov form factors to cancel 
the real soft emissions, and this resulted in a simplified expression for the gluon distri- 
bution, Eq. (I2.15p . This distribution is more "inclusive" than the original one, Eq. (12. 2p . 
which explicitly includes all soft emissions, yet it is equivalent to it for the calculation 
of the gluon distribution. In this section, we shall construct other, even more inclusive, 
versions of the CCFM evolution, which are better adapted for numerical calculations. In 
these constructions, we shall exploit the fiexibility which exists in defining the CCFM 
evolution, as associated with the various approximations involved in its derivation. Note 
that "getting more inclusive" is not the only possibility for deriving different versions of 
CCFM. In appendix B we shall derive yet another version by including non-leading effects 
related to recoils. That version could be implemented in a Monte Carlo simulation, so like 



3.1. Integral equations 

From now on we shall work with the 'unintegrated' gluon distribution, i.e., the num- 
ber of gluons with a given longitudinal momentum fraction x and a given transverse 
momentum k, which is obtained by undoing the integral over the last angular variable 
in Eq. (j2.2p (or (|2.9p ) and replacing k'^/x'^E'^. It will be also convenient to re- 

place the maximal angle ^ by a corresponding momentum variable q, via the substitution 
? = (f'/{x^E'^); as explained after Eq. (j2.2p . one has roughly ~ (the virtuality of 
the space-like photon exchanged in DIS). The integral equation satisfied by A{x,k.,q) is 



''Eventually there will be a limit coming from the soft momentum cut, so that the region does not 
extend to infinite size 




(2.22) 




(2.23) 



CASCADE. 
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easy to derive from Eq. (j2.9p . and reads 



_ /-i /■ d^p _ 
A{x,k,p) = as / dz / — :r6{p- zp)As{p,zp) 
Jx J T^P 



Ans{k,Z,p) ^ 1 



z 1 — z 



xA(^^,\k + {l-z)p\,p), (3.1) 

where we are using rescaled momenta within the integrand: p = q/{l—x) andp = q/{l — z). 
The third argument of the gluon distribution inside the integrand, i.e. p, truly means that 
the maximal angle corresponding to this distribution is the angle ^ of the emitted real 
gluon, that is, ^ = jiy^E"^) = z^p^/x^E'^, with y = (1 — z){x/z). Hence, this equation 
can be read as follows: the final t-channel gluon with energy fraction x and transverse 
momentum k (and for a maximum emission angle measured by q) is generated via the 
splitting k' k + q of a previous t-channel gluon with energy fraction x/z > x and 
transverse momentum \k'\ = \k + q\ (and for a maximum emission angle measured by p 
and z). This is the most exclusive version of the integral equation and includes all the 
hard and soft emissions. This equation is implemented in the CASCADE event generator 



In the more inclusive version we can sum over all soft emissions (in the regions Cj in 
Fig. [6|) so that the Sudakov factors disappear and we are left with the gluon distribution 
in Eq. (12.15p . This version gives rise to the following integral equation 



Aix,k,p) = asl — [ ^9ip-zp)Ansik,z,p)A(-,\k + il-z)p\,p) . (3.2) 
L z J ■np'^ \z / 



This equation is simpler to solve than (jS.ip but is still not very easy to deal with, not even 
numerically. Notice that one might as well use Eq. ()3.2p as the basis for an event generator, 
but then all soft emissions must later be included as final state radiation. In what follows, 
however, we shall be concentrating on the small-x part of the gluon distribution, where 
the small-2: values {z <C 1) are dominating. In that case one can replace all rescaled 
momenta with regular momenta, and rewrite \k + (1 — z)p\ — > |A; + g| in the argument of 
A within the integrand; also the energy ordering becomes automatic. Hence the equation 
becomes 

i^z f d'^q fx \ 

A{x,k,q) = Os / — / — K9{q- zq)Ans{k,z,q)Ai-,\k + q\,p) , (3.3) 
Jx z J Tvq^ \z J 

3.2. Geometrical representation of the real vs. virtual cancelations 

Starting with Eq. (13. 3p . we shall later derive the most inclusive, and also the simplest, 
version of the CCFM equation. To that aim it is important to understand in depth 
the structure of the virtual corrections encoded in the 'non-Sudakov' form factor Ans, 
Eq. ()2.10p . Note first that, despite its name, this form factor is essentially not different 
from a genuine Sudakov one, since it also represents the negative exponent of an area in 
the phase space, namely the areas Ai in figures [7| and El It can therefore be used to cancel 



the real emissions confined to this phase space, as was first noted in 30(]. As we review 



in appendix |A] however, this is strictly true only if the additional kinematical constraint 
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In q/k In q/k 



Figure 9: The phase space for A„s- The region B contributes with negative weight, and it is seen that 
B — B' . The shaded region A is thus the region left over, which is also the region to which the (real) 
k±-conserving emissions are conhned. To the right we show symbolically the insertion of k± -conserving 
emissions between a pair of non-k^- conserving emissions. The propagator k is unchanged by this insertion. 

k"^ > zq^, which ensures that the squared 4-momenta of the t-channel propagators are 
dominated by their transverse part [j, [i^] , is enforced within (j3.3p . 

Assuming A;^ > zq^ from now on, we can distinguish between two cases: k > q and 
k < q. When k > q, the regions Ai in figures [7] and [8] correspond to real emissions with 
q <^ k k' and A ns is guaranteed to be smaher than 1. Consider now the situation in 
Fig. [9] where k < q. As we show in appendix [Aj we have 



where the regions A (shaded region), B and B' are shown in Fig. [H The first observation 
is that the triangular regions B and B' have the same area, and hence they cancel in the 
exponent of Eq. ()3.4p . The upper diagonal line in Fig. [9] indicates the line through which 
the kinematical constraint limit k'^ = zq^ holds, while the lower diagonal line indicates 
the angle of q. Now, the emissions lying below k are, in the spirit of the approximations 
made in section [231 . kx.- conserving. When the kinematical constraint is included, these 
emissions are confined to the shaded region (region A) in figure [9l This can be understood 
as follows: due to angular ordering and the fact that we are looking at A;_L-conserving 
emissions, all subsequent radiation must lie in the region B' or in A. However, if we had 
a real emission in region i?', then because that emission is A;_L-conserving, the i-channel 
propagator emitting this gluon would have transverse momentum approximately equal to 
k and it would necessarily have bigger energy than the real gluon. Therefore we see that 
it must be located under the upper diagonal line in Fig. [9l This, however, would violate 
the kinematical constraint and such an emission is therefore not possible. Thus we are 
left with the fact that all real fcj_-conserving emissions are confined to region A. Therefore 
the inclusive summation over all the real emissions, which are inserted in between two 
non-/c_L -conserving emissions, leads to a factor 



= exp(-a,(yl + S'-S)) 



(3.4) 



= exp(asA). 



(3.5) 
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To conclude (recall that B' = B) 



A^j- A„, = exp(a,(A- A + S-S)) = 1 (3.6) 

which shows that the "non-Sudakov" factor cancels the /c i-conserving emissions. This was 
first noticed in ^(1], and it was later used in 38] as wellj. One is then left with a much 
simpler formula for A which can be derived from Eq. (|3.2p after including the constraint 
fc^ > zq^, and the fact that we are left only with the non-fc_L-conserving emissions. This 
last constraint, however, can be enforced in various ways, which are all consistent with 
each other within the present approximations. Therefore there is no unique equation that 
one can derive. In what follows we shall consider two different possibilities and then study 
the ensuing equations. 

3.3. Deriving the differential equations 

In Ref. [H we have the restriction to non~fc_L -conserving emissions by introducing the 
theta function 6{q^ — min(A;^, /c'^)) into the r.h.s. of the integral Eq. (j3.3p . This was 
also the prescription originally used in Ref. 3C], and the equations derived in [ll] and 30] 



are indeed equivalent. After inserting this constraint together with the 'kinematical' one 
k"^ > zq^ and removing A„s, Eq. (|3.3p becomes 

A{x., k, q) = a. 



[ ^ f^O{q- zq)e{e - zq^)e{q^ - mm{k\ k'^))A (-,k', q] 
Jx z J Tiq'^ \z J 



(3.7) 



Since q > k for all cases of physical interest (recall that q^ ~ in DIS), we further have 
q^ > k^ > zq^ > z^q^. Therefore the angular ordering is automatic and 9{q — zq) can 
be neglected. This means that the dependence on the third variable q drops out, at least 
in the l.h.s. But a similar argument holds also for the function for A{x/ z^k' .,q) under 
the integral, because we have q > k' (indeed, we are left only with emissions satisfying 
k' ^ q ^ k oi k ^ q ^ k'). Dropping then the dependence of A upon its third variable, 
we obtain 

A{x, k)=as f^— e{k^ - zq^)e{q^ - min(A:2, k'^))A (-, A;') . (3.8) 

J X z J 7rq \ z / 

The next step is to perform the integration over the azimuthal angle (p. To that aim, it 
is convenient to replace 9{k'^ — zq^) by 9{k'^ — zk'"^), which is allowed within the current 
approximation^, and then switch the integration variables from q to k' and, respectively, 
from z to x/ z (which we rename as z) . Making this replacement and doing the (f) integral 
one then gets 

fi dr r dk'"^ 

A[x, k) = as j ^ J |^,2_fc2| - ^fc'V^') K^i) Az, k') , (3.9) 



® In [3^, however, only the possibihties k' x k ^ q and k ~ q ^ k' were considered. In that case there 
is obviously no need for the kinematical constraint. 

^Indeed, we have either k' ~ q ^ k, in which case the replacement is obviously correct, or k ~ q ^ k' , 
in which case both the first and the second theta function can be replaced by 1. 
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where 



,M = l-^„ctan(i±^^^J.(.-l/4). (3.10) 

and K = min(fc^, A;'^)/max(/c^, /c'^). Differentiating w.r.t to Y we finally deduce the follow- 
ing differential equation for the unintegrated gluon distribution in the CCFM formalism 

r dA-'2 / 
dyAiY, k) = asj jj^^jJT\ ^(^) ( ^(^' - k'MiY, k') 

+e{k''^ - k^)e{Y - \n{k''^/k^))A{Y - ln(A;'Vfc^), k')] . (3.11) 



As mentioned above, (jS.lip is equivalent to the master equation in the LDC formalism 
(30I 1 (up to some trivial redefinitions: the distribution in 30] corresponds to our k'^A{Y, k) 
times the proton radius and some constants). 

Yet another way to implement the restriction to non-A;_L-conserving emissions is to 
switch the integration variable in Eq. (j3.3p from q to k' = \q + k\ and then replace 
with max(/c^, /c'^). Note that this constraint is more restrictive that the theta function 
in Eq. (|3.7p . In this case the angular integration in (j3.7p becomes trivial. Also the 
replacement 6{k'^ — zq^) 6{k'^ — zk''^) is now exact, and so is also the requirement 
q > k' (which, we recall, allows one to ignore the dependence of A{x, k, q) upon its third, 
'maximal angle', variable). We thus deduce 

dyA[Y, k) = as ^r^-^ - k'^)A{Y, k') 

J max(fc^,fc'^j \ 

+e{k''^ - k^)9{Y - ln{k'yk'^))A{Y - \n{k'^/k^), k')\ . (3.12) 



3.4- More on the relation to BFKL 

Before moving on to discuss the issues of unitarity and saturation, we would like 
discuss the effects of the kinematic constraint on the BFKL equation, and we would like 
to show that the same approximations used above in deriving the more inclusive equations 
(13. lip and ()3.12p applied on BFKL, with the kinematical constraint, leads exactly to the 
same equations. In section [5] we will present numerical results of the BFKL equation 
with the kinematic constraint included, with and without the saturation boundary to be 
described in the next section. 

The BFKL equation can be written 

A{Y, k) = as r dy [ ^eiY-y + \n{k'' /k'^))A'iY - y, k) A{y, k') (3.13) 

where we have included the kinematical constraint k"^ > zk'"^ (with z = e^~^). The 
function A' is the non-Sudakov form factor modified to include the kinematical constraint. 
The correct form for A' can be found by requiring that once more it can be used to exactly 
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Figure 10: The phase space for the modiSed BFKL non-Sudakov form factor when the kinematical con- 
straint is taken into account. Region B which was previously included is now excluded. Real k±-conserving 
emissions are confined to region A. 

compensate the A: j_ -conserving emissions which are now also modified by the kinematical 
constraint. The phase space for A' is illustrated in Fig. [TUJ The region we are looking for 
is the shaded one, region A. Region B which was allowed before is now excluded due to 
the kinematical constraint. Indeed the real A;_L-conserving emissions are confined to region 
A, since just like for CCFM, any real emission in B would create a t-channel propagator 
below the diagonal line, violating the kinematical constraint. It is then seen that (see also 

A'{Y -y,k)= e{k - k')A{Y - y, k) 
+ e{k' -k)exp{-asHkVQl)iHkVk'^) + Y-y)) (3.14) 

where A{Y — y,k) is the usual non-Sudakov, and the second factor is just exp{—asA) 
(again \i k > k' the kinematical constraint is automatic). 

Next we inclusively sum over the -conserving emissions in (j3.13p . canceling A'. If 
we use the explicit constraint q > min(A:, A;'), we get 

A{Y,k)=as r dy [ , „, 0{Y - y + \n{k^ /k'^))e{\k - fcf - min(A:^ k'^))A{Y, k'), 
Jo J ■n-\k-k'\'^ 

(3.15) 

and this leads to an equation which is exactly identical to (j3.9|) and (|3.11|) . We can also use 
the constraint q = max(/c, k') in which case we obtain (j3.12p . We should also mention that 
we could repeat the arguments in appendix [P] for BFKL with the kinematical constraint. 
In that case the region left over by the real-virtual cancellations is exactly equal to the 
region C in (|D.2|) . 

4. CCFM evolution vi^ith saturation boundary 

The CCFM evolution, so like any other linear evolution in perturbative QCD, predicts 
an unlimited growth of the gluon distribution with increasing Y, thus leading to unitarity 
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violations in the high energy Umit. This is so since the hnear evolution equations miss 
the non-linear phenomena responsible for unitarization, which are gluon saturation and 
multiple scattering. It turns out, however, that the phenomenon of gluon saturation 
merely acts as a kind of cutoff, which limits the growth of the gluon distribution, but does 
not modify the mechanism responsible for this growth. This makes it possible to mimic 
the effects of saturation by appropriately implementing this cutoff on the linear evolution 
equations, without a detailed understanding of the underlying non-linear phenomena. In 
this section we shall motivate and describe the implementation of this cutoff — actually, 
an absorptive boundary condition — , which then will be used, in Sect. [SI within numerical 
simulations of the CCFM and BFKL evolutions in the presence of unitarity corrections. 



4.1. Unitarity and Saturation momentum 

In what follows we shall explain our method for effectively implementing saturation 
on the example of the BFKL equation [3] . This is interesting since the corresponding non- 
linear generalization which obeys unitarity is also known — this is the Balitsky-Kovchegov 



(BK) equation — , and thus it can be used to test our method. The derivation of the 



BK equation has been recently pushed to next -to-leading order accuracy |4ll. |42|. |43|. |44||. 
but here we shall limit ourselves to its leading-order version, which is also the accuracy 
of the CCFM formalism. However, this LO version will be eventually extended to include 
a running coupling (both for BFKL and for CCFM), since the running coupling effects 
modify in an essential way the high energy evolution — these are the only NLO corrections 
which remain important for asymptotically high energy. 

The (leading-order) BK equation for the unintegrated gluon distribution reads 

n .... ,^ fdk'^(k'^A{Y,k')-k^A{Y,k) k'^A{Y,k)\ ,,,,, ,^^2 ,,,, 

The terms linear in A{Y, k) in this equation represent the BFKL equation, whereas the 
last, quadratic, term is responsible for gluon saturation. One can roughly think about 
this last term as describing the recombination of two gluons into one, but this picture is 
quite crude: the actual non-linear phenomena responsible for gluon saturation are much 
more complex and should be rather viewed as the blocking of new gluon emissions by 



strong color fields j29(]. Since this equation is non-linear, we should be more specific about 
the normalization of the function A(Y,k). Our conventions are such that the standard, 
'integrated' gluon distribution is computed as 

4/V2 /-Q^ H^A- r 
xg{x,Q^) = ^ / — / d''bA{Y,k,b), (4.2) 



vr^as J (27r)2 

where b denotes the 2-dimensional impact parameter in the transverse space. Note that, 
strictly speaking, the non-linear effects are non-local in b. The simple form for the non- 
linear term shown in Eq. (|4.1|) is obtained under the further assumption that the hadron 
is homogeneous in b (a 'large nucleus'). With these conventions, the gluon occupation 
number — i.e., the number of gluons of a given color per unit rapidity per unit volume in 
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transverse phase-space — is not exactly A(Y, k), but rather A{Y, k)/as (up to a numerical 
factor) . 

So long as A(Y,k) <^ 1, the non-linear terms in Eq. ()4.ip can be neglected, and 
then this equation predicts the rapid, BFKL, growth of the gluon distribution, which is 
exponential in Y, together with diffusion in transverse space (see below). When A(Y, k) ~ 
1 (corresponding to a physical occupation number 0{i/6is)), the non-linear effects become 
important and tame this growth. For a given rapidity Y , this happens at a specific value 
of the transverse momentum k = QgiY), called the saturation momentum, which grows 
rapidly with Y (see below). One can show that at low momenta k < QgiY), the occupation 
number essentially saturates, in the sense that is shows only a very slow increasqH with 
Y: AiY, k) ~ ln[Qs{Y)/k]. Hence, for a given Y, the gluon distribution A{Y, k) produced 
by the BK equation looks like a front, which interpolates between the dilute {A{Y, k) <^ 
1) BFKL taiO at high transverse momenta k S> Qs{Y) and a saturation region at low 
transverse momenta k < Qs(Y), where A{Y,k) ~ ln[Qs(y)//c]. The transition between 
these two region occurs around k = Qs{Y), where A{Y,k) ~ 1. With increasing Y, this 
transition value Qs{Y) is rapidly increasing, i.e., the front moves up to higher values of k. 

What is remarkable about this dynamics is that the progression of the front with 
increasing energy and also its shape at high k are fully determined by the BFKL evolution 
of the dilute tail at k ^ Qs(Y), and thus can be inferred from the linear, BFKL, equation 
alone. One says that the saturation front is "pulled by its tail". This property is central 
to our analysis: it implies that some essential features of the dynamics in the presence 
of saturation, like the energy dependence of the saturation momentum, can be studied 
without a detailed knowledge of the non-linear effects responsible for saturation. Hence, 
a similar study can be performed on the basis of other equations, so like CCFM, whose 
non-linear generalizations are not known. 

The pulled-front property property is highly non-trivial — there are many examples of 
non-linear equations which develop a pushed front, i.e., a front whose progression is driven 
by the growth and accumulation of 'matter' behind the front 45|] — and so far it has not 
been rigorously demonstrated for the general case in QCD. In the case of a fixed coupling 
and for sufficiently high energy, this property follows from the identification [8] between 
the asymptotic form of the BK equation at high energjifl and the FKPP equation (from 
Fisher Kolmogorov, Petrovsky, and Piscounov) of statistical physics. (The FKPP equation 
describes a 'reaction-diffusion process' in the mean field approximation corresponding to 
very large occupation numbers at saturation; see e.g. the review paper [i^.) However, 
this identification does not extend to a running coupling, and it was in fact shown 
that the high-energy evolution with running coupling is not in the same universality class 
as the reaction-diffusion. Yet, the pulled-front property appears to hold for that case 



^The dependence of Qs{Y) upon Y is such that the physical occupation number at fc < Qs{Y) grows 
hnearly with Y; that is, what saturates is the rate for gluon emission [29| . 

*In particular, for extremely high k 2^ Qs{Y), the BFKL distribution approaches the bremsstrahlung 
spectrum, A{Y,k) ~ whereas at moderate k this is modified by the BFKL 'anomalous dimension'. 

^This is obtained via the gradient expansion of the non-locality in Eq. (|4.ip to second order in d/dp, 
with p = Infc^ ('diffusion approximation'). 
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too, and also for relatively small rapidities {i.e., for the early stages of the evolution), as 
most convincingly demonstrated so far by the numerical simulations in Ref. [H , where the 
solutions to the BK equation with running coupling have been systematically compared 
to those of the BFKL equation supplemented with a saturation boundary condition. 

In order to gain some analytic insight into the role and the form of this boundary 
condition, it is useful to briefly review the computation of the saturation momentum from 
the BFKL equation 0, 0] — in the fixed coupling case, for simplicity (the corresponding 
developments for the case of a running coupling can be found in Refs. 0, [IB])- We start 
with the Mellin representation of the BFKL solution, that is 

A{Y,k) = y ^ e"»^(^)^-(i-^)M(7), (4.3) 
c 

where p = ln{k'^ /Qq) with Qq an arbitrary reference scale, ^(7) is the initial condition 
at y = 0, and xil) is the BFKL characteristic function, i.e. the eigenvalue of the BFKL 
kernel in Mellin space: 

x(7) = 2V(l)-V'(7)-^(l-7), ^(7) = dlnr(7)/d7. (4.4) 

The integration contour C runs parallel to the imaginary axis with < Re(7) < 1. For 
real values of 7 in between and 1 (the relevant range for computing the saddle point; 
see below), the function xil) is displayed in Fig. [TTJ 

Eq. (I4.3jl is expected to be correct so long as A{Y,k) <^ 1, i.e. for high enough 
momenta k ^ Qs{Y). Here we shall assume that this expression can be also used to 
approach the saturation line 'from the above' (i.e. from momenta k larger than Qs), and 
hence to approximately determine the latter from the condition that A{Y, A;) ~ 1 when 
k = Qs{Y). To that aim, we shall also rely on the saddle point approximation, which is 
appropriate for high enough Y — namely, such that agY ^ 1. By combining the saddle 
point condition 

asYx'{-fs) = -Ps, (4.5) 
with the condition that A{Y, k) 1 along the saturation line: 

asYxils) - (1 - ls)Ps = 0, (4.6) 
one obtains 7^ and the saturation line PsiY) = ln[Q'l{Y)/kQ] as follows 
X'ils) _ 1 



Xils) 1-7 



7s ~ 0.372, 



s 



Ps{Y) ~ XsY with A, = « 4.883a, . (4.7) 

1 -7s 

A geometrical interpretation of the above equation for 7^ is shown in Fig. [TTJ 

For p larger than ps, but not much larger, one can estimate the gluon distribution by 
expanding the integrand in Eq. (j4.3p around 7, : writing 7 = 7s — ii^, and expanding to 
second order in v, one finds 

asYxil) - (1 - 7)P ^ -(1 - 7s){p - Ps) - ^KP - Ps) - D^Yv^ (4.8) 
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Figure 11: The BFKL eigenvalue xil) ^-nd the graphical solution to the saturation problem. The value 
7s = 0.327 corresponds to the saturation saddle point and the value x(7s)/7s = 4.88 determines the 
asymptotic energy dependence of the saturation momentum (for comparison, the "hard pomeron" saddle 
point 7p = 1/2 and its intercept = 41n2 = 2.77 are shown). 



where Dg = asx" {is) /'^ ~ 24.260^. This expansion is vahd so long as 1 < p — ps DgY. 
After also performing the Gaussian integration over u, one finds [g] 



with Q^{Y) = Q^e^"^ . Eq. (j4.9p exhibits a power-like spectrum in k with anomalous 
dimension (i.e., deviation from the bremsstrahlung spectrum) 7^, which is further modified 
by BFKL diffusion (the last factor, which is a Gaussian in ln(fe^/Q^) = p — Ps)- 

Although obtained solely from the linear, BFKL, equation, the above results for the 
behaviour near saturation are essentially correct: they coincide with the respective pre- 
dictions of the BK equation for asymptotically high energy. It turns out that one can 
do even better: still without resorting on the actual non-linear Eq. (j4.ip . one can also 
determine the sub-asymptotic behaviour of Ql(Y) at large Y, and refine the approach 
of A{Y, k) towards saturation. The main observation, which lies also at the heart of the 
subsequent developments in this work, is that the saturation region acts like an absorptive 
boundary [3], which not only tames the growth of the gluon distribution, but also prevents 
the BFKL diffusion towards lower momenta k < Qs{Y). 

After also implementing this absorptive boundary, in a way to be shortly specified, 
the previous formula are modified as follows [7|]: the gluon distribution becomes 

A{Y,p) ~ (p-p, + 5)e-(i-^=)(^-''=) expj-^^^^j, (4.10) 
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where (5 ~ 1 is some unknown constant and PsiY) is now under control up to next -to- 
leading order in the asymptotic expansion at high-energy: 



Ps{Y) ~ XsY - -^Iny. 



(4.11) 



's 



These results coincide indeed with the corresponding predictions of the non-linear, BK, 
equation, as deduced via the correspondence with the FKPP equation 0]. Notice in 
particular the property of geometric scaling 0, 0]: within the region p — ps <^ y/lDgY , 
whose width is increasing with Y , the Gaussian factor can be ignored in Eq. (|4.1U|) . and 
then this expression reduces to a function of p — Ps{Y) alone, i.e. of the scaling variable 
Q"^ /Q1{Y). This property provides a natural explanation for an important regularity 
observed in the HERA data for DIS at small-x 1161 1. 




The success of the absorptive boundary method for the BFKL evolution makes it 
compelling to try and use it as a systematic method for enforcing saturation within an 
arbitrary linear evolution, so like CCFM. Let us now explain in detail our practical im- 
plementation of the absorptive boundary, as it will be used within numerical simulations. 
We first introduce a line of constant gluon occupancy p = Pc{Y) via the condition 



where the number c is smaller than one, but not much smaller. The saturation line PsiY) 
would correspond to c ~ 1, so clearly PsiY) is smaller than p^iY), but relatively close 
to it. (For the BK evolution with fixed coupling, these two lines would be parallel lines, 
with slope As, separated from each other by PciY) — Ps{Y) ~ ln(l/c).) For p < Pc{Y) and 
sufficiently high energy, the solution Abfkl{Y, p) to the BFKL equation would become 
larger than one — in fact, arbitrarily large. If this equation is to be solved numerically, 
one may think about enforcing saturation by hand, in the following way: at each step in Y, 
one first identifies the corresponding point pdY) from the condition ()4.12p . and then one 
requires A{Y, p) to remain finite and of 0(1) for any p sufficiently far below PciY) — say, 
for p < PciY) — A with A ~ ln(l/c). When decreasing p below pdY), the solution A{Y, p) 
will typically start by rising, then reach a maximum of 0{\) and eventually decrease to 
zero. We shall conventionally identify the saturation scale Ps{Y) with the position of this 
maximum. In this procedure, the numbers c and A are to be viewed as free parameters, 
which are however correlated with each other, since A ~ ln(l/c). 

By construction, the value of AiY, p) behind the saturation front is not under control, 
as this is fixed by hand at a constant value. In practice we shall choose this constant 
value to be zero: A{Y-:P) = for p < Pc{Y) — A (but other values of 0(1) will be also 
used, to test the sensitivity of the method to this particular choice). Thus, clearly, our 
procedure cannot be used for those physical problems which are sensitive to the details 
of the saturation region, like deep inelastic scattering at low < Ql(Y), or particle 
production at low transverse momenta. On the other hand, this procedure accurately 
describes the dynamics of the front, in that it provides the same results as the non-linear 
BK equation for the energy dependence of the saturation momentum and for the gluon 
distribution AiY^ k) at momenta k > QgiY)^ for both fixed and running coupling, and for 
all values of Y . 




A{Y,p = pdY)) 



= c 



(4.12) 
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Let us illustrate the efficiency of this method with a few numerical results for the case 
of a running coupling. We include running coupling by pulling the ots factor inside the k'- 
integral in Eq. (14. ip and using the one-loop expression for the running coupling with scale 
= max(fc^, k''^) and Aqcd = 200 MeV. This simple prescription is in agreement with the 
recently constructed running-coupling version of the BK equation [4l|, [Z^] . To avoid the 
infrared divergence of the coupling at = ^qqd, we shall replace as{Q'^) — > ctsiQ^ + 
for some parameter /i. Our default choice will be /i^ = 0.5 GeV^, but we shall study the 
sensitivity of our results to variations in ^. Our initial condition A(Y = 0, k) is given 
by the bremsstrahlung spectrum for k > 1 GeV (with maximal height A = 0.5) and it 
vanishes for < 1 GeV. 

The pure BFKL evolution with running coupling is known to be infrared unstable: 
the rise of the gluon distribution is much faster at small values of k (where the coupling is 
larger) and, besides, this rapid accumulation of gluons in the infrared is also feeding the 
growth at higher k. Therefore the linear evolution behaves quite differently compared to 
the non-linear one, even at high k. To illustrate that, we have also included in Fig. [12] 
the results of the strict BFKL evolution, which indeed show a much faster progression 
towards high k as compared to the BK equation. On the other hand, one observes a 
perfect matching between the saturation fronts provided by BK and, respectively, BFKL 
with saturation boundary. This shows that the front remains of the pulled type even 
with running coupling and at the same time demonstrates the success of our method for 
effectively implementing saturation. We also note that in the corresponding figures in [2], 
the values c = 0.1 and A = 5.0 were used instead of c = 0.4 and A = 2.0, thus showing 
that the results do not depend on the specific values of these parameters (as long as they 
are correlated as A ~ ln(l/c)). 

In particular, the infrared problem is cured by saturation: the saturation scale effec- 
tively acts as an infrared cutoff, which becomes 'hard' {Q'l{Y) ^ ^qcd) sufficiently 
high energy. To better illustrate this, we exhibit in Fig. [13] results obtained for different 
values of the IR cutoff //^ inserted in the running coupling. Unlike the pure BFKL results 
(left figure), which are extremely sensitive to a change in /i, the results corresponding to 
the saturation boundary condition (right figure) show no sensitivity whatsoever. 

From the curves in Fig. [T2l it is also possible to extract the ^-dependence of the 
saturation momentum Ps{y) for running coupling. We find that the squared-root law 
XrVy predicted by the theory 0,0, 15 1 for asymptotically high energies provides a 



good fit to our numerical results for Y > 10, with a fitted value Ar- — 2.9 which agrees 
reasonably well with the (asymptotic) theoretical expectatior0 A^. — 3.2. 

Now, from the phenomenological point of view, we are more interested in values of Y 
which are not that large, say y < 14 (corresponding to x > 10^^), as relevant for forward 
jet production at LHC. With that in mind, we also show in Fig. [12] (right) the results for 
lower values of Y, between 6 and 12 units; one can thus see that the absorptive boundary 



For asymptotically large Y, the running-coupling BFKL evolution yields I, [3] : PsiY) ~ V2\oboY 
where Aq — 4.88 is the same number as in Eq. (|4.11|l and feo = 12Nc/ (llNc — 2Nf) is the coefficient in 
the one-loop running coupling: q(Q ) — ho/\n{Q /Aqcd). In our simulations, we use Nf — 0, hence we 
expect Ar = V2Ao6o — 3.26, which is indeed consistent with the fit to the curves in Fig. 1121 
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k (GeV) k (GeV) 

Figure 12: The running coupling results for: BK (solid curves), BFKL with absorptive boundary (long 
dashed curves) and pure BFKL (short dashed curves) for (left) Y = 10,20,30 and 40, and (right) Y — 
6, 8, 10 and 12. For the absorptive boundary we used c — 0.4 and A = 2.0. 
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Figure 13: Testing the sensitivity of the BFKL evolution to the IR cutoff ^. Left: The pure BFKL evolution 
with ^? = 0.5 GeV^ (solid lines) and = 2.0 GeV^ (dashed lines), for Y = 10, 20 and 30. Right: The 
BFKL evolution with absorptive boundary for — 0.5 GeV^ (solid lines), fi'^ — 2.0 GeV^ (dashed lines) 
and for the same values for Y as before. 

method works equally well also for such lower rapidities. 

4-2. CCFM evolution in presence of saturation: analytic results 

Within numerical simulations, it is straightforward to implement the saturation bound- 
ary condition in any of the previous versions of the CCFM equation: the complete, integral, 
equation (13. ip . or one of its simpler, differential, equations deduced in Sect. 13.31 Given 
the approximations made in deriving the latter, it is clear that the corresponding results 
will not be exactly the same. In particular, the saturation momentum deduced from the 
full CCFM equation (j3.ip will also depend upon the maximal angle variable ^ introduced 
by the kinematics of the external scattering; in the case of DIS, this means that Qa will 
depend upon the virtuality scale at which one measures the structure functions. It 
would be very interesting to study this dependence, and also the energy dependence of 
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Qs, by numerically solving the integral equation (j3.ip in the presence of saturation. This 
could be done via Monte Carlo simulations (e.g., by correspondingly extending the CAS- 
CADE event generator which would also allow one to study the effects of saturation 
on the structure of the final state. Such numerical calculations would be however tedious 
in practice, and it would be very difficult to reach high values of Y in this way. Our 
interest here is mostly conceptual: we would like to understand how saturation modifies 
the CCFM evolution, and how the respective predictions compare to those based on the 
BFKL evolution (with saturation once again). To that aim, we shall focus on the differ- 
ential versions of the CCFM equation, for which we shall present some analytic estimates 
in this section, and then a systematic numerical study in Sect. [5l 

The analytic study is based on the Mellin transform, as already introduced in Eq. (14. Sp . 
In what follows we will generically denote the kernel eigenfunction (the quantity denoted 
as Oisxil) ill (14.31) ) as w(q;s, 7), keeping the notation ^(7) only for the BFKL eigenfunction 
(14. 4p . In the case of Eq. ()3.12p . the function uj can be analytically computed, while for 
Eq. (j3.1ip it must be constructed numerically. Namely, for Eq. (j3.12p . we get 



u = as[- + ) , (4.13) 

7 1 — 7 + 



which is easily solved to give 




/ — \ I / — \ 1 — 



For as ^ this expression reduces to 

io = as(- + (4.15) 
V7 1-7/ 

where I/7 is recognized as the collinear piece (the one which dominates when k ^> k'), 
while 1/(1 — 7) is the anti-collinear one {k' ^ k). The high-energy saddle point coincides 
with the BFKL one, 7 = 1/2, but the intercept is given by 4, instead of 2.77 for BFKL. 
The saturation saddle point, 7^, is found according to Eq. (j4.7p which now becomes 

-'(7.) = . (4.16) 

1 - 7s 

This gives 7^ = 1/3 ~ 0.33, which is quite close to the respective BFKL value 7^ ~ 0.37. 
The saturation exponent controlling the growth of the saturation momentum is again 
obtained from (j4.7p which in this case gives 6.75, and thus is larger than the respective 
BFKL value 4.88. 

Although the eigenfunction ()4.15p gives a too high intercept and a too large speed of 
the saturation front, it is well known that it approximates the general shape for the BFKL 
eigenfunction, x(7)) very well. (This explains why the respective saturation saddle points 
are so close to each other.) In fact if one just subtracts from (j4.15p the difference of the 
two intercepts, 4 — 41n2, then an almost perfect approximation of xil) is obtained. Note 
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Figure 14: Left: The Mellin space eigenfunction of the asymptotic approximation of Eq. (|3.1ip (soUd line) 
compared to the BFKL eigenfunction xil) (dashed line). Right: The same as in left but with a constant 
0.45 subtracted from the solid line. 



that a constant shift in the eigenfunction corresponds in momentum space to the inclusion 
of a local term oc A{Y, k) in the right-hand side of the differential equation. 

For Eq. (j3.1ip . we write u = ctsxil,'^), where x can be numerically constructed. In 
the weak coupling limit (q^ 0) we have 

u; = asxh) = dt ^(t'^ + f-^'-^^). (4.17) 

JO ^ — 1^ 

To find the behaviour at the endpoints, 7 = and 7 = 1, we concentrate on the region 
t ^0 and use that h{t) = 1 for t < 1/4. This yields 

1/4 1 

l-V ^ 

= ^ ; ' F{1 - 7, 1, 2 - 7, 1/4) + ^^F(7, 1,1 + 7, 1/4), (4.18) 
1-7 7 

where F is the hypergeometric function. Since at the endpoints we have F ^ 1 for each 
term, we find the same endpoint behaviour as ()4.15p . and hence as in BFKL. For the 
complete interval we plot x, together with x, in Fig. [TH We see that the minimum of x, 
i.e. the intercept, is a bit above the BFKL one — this is numerically found as 3.23 — , and 
therefore so is also the respective saturation exponent, found as — 5.56. On the other 
hand, the shape of the eigenfunction is again very similar to the BFKL one: in the right 
figure [T^ we have subtracted the difference between the two intercepts, i.e. 0.45, from x, 
and then the two curves almost coincide. So, no surprisingly, the saturation anomalous 
dimension 7^ ~ 0.35 is close to the respective BFKL value. As previously noted, such 
a constant subtraction is tantamount to a contribution proportional to A(Y, k) in the 
differential equation, which in turn can be associated with additional virtual corrections. 
In appendix [P] we will see that such virtual terms naturally appears in Eq. (j3.1ip when 
the real-virtual cancellations are treated with a better accuracy. 

Before leaving this section, we demonstrate the behaviour of the eigenfunctions for any 
as- For finite a^, the anti-collinear pole is screened by uj appearing in the r.h.s. of (|4.13p 
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Figure 15: The eigenfunctions u)/as for equations (|3.11|) (solid line) and (|3.12p (dashed line), for 4 values of 
the coupling (from up to down) : as = 0.001 , 0.01, 0.1 , and 0.4. The very small values of ds are included 
to illustrate the emergence of the anti-collinear pole at y = 1 in the limit Us 0. 

and the corresponding equation for (j3.11|) . The pole appearing in (j4.15|) is really valid 
only when = 0. In Fig. [15] we plot the eigenfunctions, uj/ctg, for (jS.lip and ()3.12p for 
various values of dig- The figure clearly demonstrates the emergence of the anti-collinear 
pole as as — > 0. 

5. Numerical results 

In this section, we shall present numerical solutions for the two differential versions 
of the CCFM evolution derived in section [3.31 namely Eqs. ()3.1ip and (I3.12p . and also 
for the BFKL equation including the kinematical constraint, cf. Eqs. (|3.13p and (j3.14p . 
We shall mostly focus on the unitarity-preserving versions of these equations, as obtained 
after implementing the saturation boundary described in Sect. I4.1i But before doing that, 
it is also instructive to compare the respective linear evolutions, in order to illustrate 
the significant differences between them (in agreement with the discussion in Sect. 14. ip . 
and thus emphasize the ambiguity in our current understanding of high-energy evolution 
within perturbative QCD. Remarkably, it will turn out that this ambiguity is drastically 
reduced after including the saturation effects in the form of the absorptive boundary. 
Namely, the main conclusion that will emerge from our analysis is that all the equations 
under consideration — the CCFM equations (j3.1ip and (j3.12p . and the BFKL equation 
with and without the kinematical constraint — , which differ significantly from each other 
in the linear regime, lead nevertheless to very similar predictions for the energy dependence 
of the saturation momentum after adding the saturation boundary condition and with a 
running coupling. Note that the inclusion of the running coupling is essential in that sense: 
with a fixed coupling, important differences persist even after including saturation, as it 
should be expected from the discussion in Sect. 14.11 
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Figure 16: The solutions to 1)3. lip (solid lines), and the BFKL equation including the kinematical constraint 
(|5.1|) (dashed lines) for Y — 10, 20 and 30. Left: fixed coupling as = 0.2. Right: running coupling. 

For the purposes of the numerics, it is convenient to rewrite Eq. (j3.13p in differential 
form and perform the azimuthal integration; this yields: 



In Fig. [16] we show the solutions to the linear equations (|3.11|) and (j5.ip . As expected 
from the discussion in Sect. 14. ll one sees that the CCFM solution is somewhat faster. The 
difference is more pronounced for a running coupling, due to the infrared instability of 
the linear evolution. But even with a running coupling, the present equations (15. ip and 
(13. lip , which include the kinematical constraint, predict a growth for the gluon distribution 
which is strongly reduced as compared to the strict BFKL equation (without kinematical 
constraint). This was shown in Ref. [H where the solutions of ()3.1ip have been compared 
to the solutions of the pure BFKL equation. 

Let us now perform a similar comparison after adding the saturation boundary. (Un- 
less otherwise stated, we use the values c = 0.1 and A = 5.0 for all calculations including 
the saturation boundary. A different set of values will be used later on, to illustrate the 
robustness of our results.) For the same Y values as in Fig. [TBI we compare in Fig. [T7] 
the respective predictions of Eqs. ()3.1ip and ()5.ip for both fixed and running coupling. To 
make the comparison between the solutions easier, we have scaled the solutions to (j5.ip 
so that the solutions more or less coincide at Y = 10. For the fixed coupling case (left 
plot) the scaling factor is 3.9 while for the running coupling case (right plot) it is 3.6. The 
difference between the fixed and the running coupling cases is striking. For fixed coupling, 
the saturation front generated by ()3.1ip is seen to progress significantly faster than the 
one generated by (j5.ip . For the running coupling case on the other hand, it is clear that 
the respective fronts progress with similar speeds. Thus in this case the two equations 
have similar predictions for the energy dependence of the saturation momentum, so that 
the entire difference between them can (almost) be removed by either a rescaling of the 



dyAiY, k) = a J dk' 
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Figure 17: The solutions to equations 1)3. lip (solid lines) and (15. ip (dashed lines) including the saturation 
boundary, for Y — 10, 20, 30. The solutions to (|5.ip have been scaled so that the solutions match at 
Y = 10. Left: fixed coupling as — 0.2. Right: running coupling. 
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Figure 18: Solutions to (|3.1ip . (|5.1I) and BFKL including the saturation boundary for asymptotic Y 
between 50 and 140 units, and for a running coupling. Left: (|3.11l) (dashed lines) and (15.11) (solid lines). 
The solution to (|5.ip has been scaled by a factor 6.5. Right: BFKL (dashed lines) and (|5.ip (solid lines). 
The solution to ()5.ip has been scaled by a factor 4. 



saturation momentum Qs, or as in the figure, a rescaling of A. This is clearly not the case 
for a fixed coupling. Note also that the difference in the fixed coupling case is visibly large 
even though the scale of the figure is much larger than corresponding running coupling 
plot. 

In Fig. [18] (left) we show the respective results for the running coupling case for 
much higher values of Y. Also, in the right figure there, we compare the absorptive 
boundary results for Eq. (|5.1|) and, respectively, the standard BFKL equation without the 
kinematical constraint. One clearly sees now that, as anticipated, all the equations under 
consideration give approximately the same speed for the saturation front, at least up to 
Y = 140. To make this similarity more manifest, we have rescaled the gluon distributions 
by appropriate factors, in such a way to superpose the different sets of fronts. These 
factors account for the difference in the overall normalization of the saturation momenta 
provided by the different equations. 
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Figure 19: The solutions to equations (|3.1ip (solid lines) and p.l2|) (dashed lines) in the linear case with 
a running coupling for Y = 10, 20 and 30. 




k (GeV) k (GeV) 

Figure 20: The solutions to equations (jS-lip (solid lines) and (|3.12p (dashed lines) including the sat- 
uration boundary and for a running coupling for Left: Very high Y between 40 and 100. Right: The 
phenomenologically relevant values ofY= 10, 12, 14. 

We have also checked that the energy dependence of the saturation momentum is 
consistent with the form expected at running couphng, i.e., = QqC^''^. Fitting this 
form to our results we find A,- ~ 3.1 for Eq. (|3.1ip . and ~ 2.9 for Eq. (jS.ip . The 
normalization factor Qq is found to be a factor 1.6 higher for (13. lip then for ()5.ip . For 
BFKL without the kinematical constraint we again find A^ ~ 2.9. 

Figure [19] shows the comparison between the two versions of the CCFM evolution, 
Eqs. (|3.1ip and (j3.12p . in the linear case (no saturation boundary) and for a running 
coupling. We see that (|3.12p gives a faster growth, but the shapes of the curves are 
similar; this was to be expected in view of the discussion in Sect. 14. 2i The corresponding 
results after adding the saturation boundary are shown in Fig. [20l It is then find that the 
speeds of the saturation fronts are quite similar, with Eq. ()3.12p giving a slightly larger 
speed, meaning a higher value for the saturation exponent. Up to the Y values shown in 
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k (GeV) k (GeV) 

Figure 21: Comparing the solutions to Eq. (|3.11|) (with saturation boundary and running coupling) using 
two set of values for the saturation parameters: c — 0.1, A — 5.0 (solid lines), c = 0.4, A = 2.0 (dashed 
lines). Left: Y = 8, 10, 12, 14. Right: Y between 30 and 100. 



the figure, we find the value Ar ~ 3.3 for ()3.12p compared to the value Ar ~ 3.1 for ()3.1ip 
mentioned before. The difference in the normalization factor Qq is a factor 1.4, with (j3.12p 
giving the larger value. 

So far all these results have been obtained for the fixed set of values c = 0.1 and 
A = 5.0. Theoretically, the different values of c and A should only affect the absolute 
normalization of the gluon distribution and of Qs, while leaving unchanged their shape and, 
respectively, energy dependence. To check these expectations and illustrate the robustness 
of our results we now repeat some of the calculations with the set c = 0.4 and A = 2.0. 
(Recall that these two parameters are correlated as A ~ Inl/c.) 

In Fig. [21] we compare the results obtained with the two sets of values for the case of 
Eq. (j3.1ip . For the phenomenologically relevant values between Y = 8 and 14, it appears 
that the curves are moving with slightly different speeds, and it seems that there is a 
slight difference in the slope of the gluon distribution which cannot be removed by a pure 
rescaling. The differences are, however, tiny and are not significant within the numerical 
certainty. For a much larger interval Y where the differences should become more visible, 
we in fact see that the shapes of the curves and their speed are very similar to each other. 
The difference in normalization is around a factor 1.5, with the solutions with c = 0.4 
and A = 2.0 lying below the default ones. Furthermore, in Fig. [22] we show a comparison 
between Eqs. (j3.1ip and (|5.ip using this new set of parameters. We have here again scaled 
up the solutions to (j5.ip . but this time by a factor of 5.5, compared to 6.5 in Fig. [18] Once 
again, the two sets of fronts appear to have similar speeds and shapes, with the CCFM 
solution progressing slightly faster than the BFKL one. We have checked that the speeds 
of the fronts come out the same as in the default case of Fig. [18] 

Let us finally consider different versions of the saturation boundary, to see whether 
this might affect our results. So far we have always applied a totally absorptive boundary 
to mimic the non-linear physics of saturation. From a theoretical viewpoint, it should be 
no fundamental difference between enforcing A to vanish behind the saturation front, or 
fixing it to some non-zero constant value of order one: various choices for this value must 
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Figure 22: Comparing the solutions of Eqs. (|3.1ip and (|5.1[) for c — 0.4, A = 2.0, and large Y between 30 
and 100 units. The solutions to (|5.ip have been scaled up by a factor 5.5. 
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Figure 23: The solutions to Eq. (|3.11[) with different boundary conditions for a running (is, and for 
Y = 10, 20, 30, 40 and 50. Left: A = 0.5 (solid lines) compared to the default A = (dashed lines). Right: 
A = 1.5 (solid lines) compared to the default ^ = (dashed lines). 

lead to the same front dynamics at asymptotically high energies. One may however wonder 
whether differences can be important in practice, for non asymptotic energies. We have 
therefore tried two other boundary conditions as well, namely A = 0.5 and respectively 
A = 1.5. The corresponding results are shown in Fig. [23l where they are also compared to 
the results obtained with the absorptive boundary condition ^ = 0. The differences are 
find to be rather minor and correspond more or less to rescaling the gluon distribution. 
In particular, the energy dependence of Qs and the shape of A come out the same in all 
cases and for values of Y, including the relatively small ones. 
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Figure 24: The representation of the situation where the real gluon (gluon b) takes a large fraction of 
the transverse momentum of the virtual propagator emitting it ( Qa)- The next propagator Qt satishes 
Qb < Qb- In this case the correct non-Sudakov is obtained by integrating over region A with opposite sign. 
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A. On the correct form of Ans 

In this appendix we would like to make clear the correct form of the non-Sudakov 
form factor Ans, as it is usually written in a way which is not entirely correct. In the 
literature one usually finds the formula 

- exp(-.,,n(5),„(-|L)). (A., 



where 



1 if Qk/Pk > 1 
Qk/Pk if Zk < Qk/Pk < 1 
Zk if Qk/Pk < Zk 



In this formula, the integrand is always negative so that A„s always gives a suppression. 
Although this seems reasonable, it is, however, not correct. 

When we identified the areas in Fig. [7] we assumed Qk > z^qu- Consider now the 
set of emissions in Fig. [211 Here we also show the maximum allowed angle ^. Note that 
the region B would have been included in BFKL but is excluded in CCFM. The smaller 
regions c, d and e are included in CCFM. The interesting region here is the one marked 
by A. It is contained within the region Ch and therefore also in the Sudakov associated 
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In 1/z 




Figure 25: Closer inspection of the phase region over which the non-Sudakov is integrated. The correct 
result is obtained by subtracting the regions B (in figure (a) ) and A (in Bgure (b)). The results are given 
in the text. 



with the gluon b. When b is emitted from Qa, it takes almost all transverse momentum 
so that Qa ^ Qb ^ Qb- The region A is thus bounded from below by Qf, as shown in the 
figure. Here indeed ZbQb > Qb, and if we set A„e = 1 here so that it does not contribute 
then we have only the contribution from the Sudakov form factor in this region. However, 
there should actually be no net contribution from the form factors in this region. To see 
this one should remember the result in Fig. [3l when we multiply Sne and Seik, what is left 
over is a region bounded from above by the virtual propagators Qk- 

There should therefore be no net contribution from the form factors in region A. 
Thus in this region the Sudakov and the non-Sudakov must cancel each other. If this is 
to happen we have to set 

Ane = exp i+aA) (A.2) 

here and not Am = 1- 

Actually let us from geometrical considerations derive the formula for A^s more care- 
fully. In Fig. [25] we show the phase space over which the non-Sudakov is integrated. 
In figure (a) we have (denoting the virtual gluon momentum by k and the real gluon 
momentum by q) 

Ans = exp(-a,(5 - A)) (A.3) 
where the areas A and B are marked in the figure. We see that 

B-A = ln'--ln^- = ln-ln^. (A.4) 

zq q z zq^ 

For figure b we instead have 

= exp(-a,(-^)), (A.5) 
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and 

^ = lniln-f- -ln2- = -In-ln^, (A.6) 
z z z zq'^ 

while for figure (c) we have 

A„, = exp(-a,^), (A.7) 
and ^ ^ 

^ = ln-ln^+ln2- =ln-ln-^. (A.8) 

z z z zq^ 

Thus A„s is indeed given by the formula in (|2.2ip . As remarked before we see that A^^ 
always gives a suppression if k"^ > zq'^, that is if the kinematical constraint is assumed to 
hold. 



B. Non leading effects in the CCFM equation 

When deriving the formulas for A„s and A^ in (I2.2ip and (12.230 from the original form 
factors in (j2.3p and (|2.4p . recoils were not taken into effect. In this section we will propose 
modifications to ()2.2ip and ()2.23p in order to more properly take into account recoil effects. 
Although such effects are formally suppressed, experience tells us that they can be quite 



important for phenomenology. For example in reference 46|], Ans was modified beyond 
the leading order in a simple way, and it turns out that the modification is significant for 
phenomenology. Besides, in the Monte Carlo implementation in CASCADE full energy- 
momentum conservation is already taken into account and one has then no reason not to 
modify the virtual form factors accordingly. 

For the sake of demonstration we will distinguish between the momenta of the t- 
channel gluons in the angular and, respectively, energy ordered cascades, that we shall 
denote as Ki and, respectively, Qi. (Recall that Ki and Qi differ by soft gluon effects 
which have been neglected in the main text.) In Fig. [26] we again show a gluon chain, 
but this time including the recoils that were so far neglected. The shaded areas here show 
the difference in phase space between the original expressions in (|2.3p and (|2.4p . and Ans 
and A^. The horizontal dashed lines in the figure represent Qi while the Ki are marked 
by crosses as before. So actually also the region between the dashed horizontal lines and 
the solid lines represent the differences although we have not shaded them. The problem 
is that it is difficult to take into account this difference since we cannot construct the 
exact Qi until the whole chain has been generated, as the final energy ordering is a priori 
not known. We will, however, see how we can take into account the shaded regions by 
modifying A^^ and A^. 

What we see is that one generally overestimates the phase region included in A^s 
and As- This means that one overestimates the suppression coming from these factors. 
Actually in some of the shaded regions in the figure both Ans and overestimate the 
phase region so we get double extra suppression. Let us first concentrate on the soft 
emissions. For the soft emission pi, the y integral is integrated up to Xi-i. Now this will 
always cause an overestimate since the subsequent emissions are emitted from Ki with 
energy Xi < and when eventually the next hard gluon is emitted it will have energy 
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Figure 26: A CCFM chain in the phase space, including also recoils in the emissions. The shaded regions 
indicate the difference in phase space between the original expression in (|2.3|) and (|2.4p . and the form 
factors Ans and As . 



less than Xi itself. Thus we see that it would be better to have the energy integral in 
for the soft emission integrated up to Xi, instead of Xj-i. Ideally we would like to have 
had the integral up to the energy of the next hard gluon, but of course we do not know 
what that energy will be beforehand. Thus we propose that 

Ar^*(z) ^ exp (-a. ^ -\ (B-l) 




for soft emissions. Again writing this in terms of momenta we get 




which should be compared to (12.230 (again e' = qo/p)- We see that the difference in the y 
integral is 

y Je' y Zi 

where of course for soft emissions ~ 1. 

However, now we can also see that there is a region between and Xi which is 
missed since the contribution from A^"-^* has been removed from this region, as illustrated 
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overestimated 




Figure 27: The shaded regions marked "missed" are regions which contribute to the original form factors 
but not to the new form factors, after the Sudakov for the soft emissions has been moBfied. Regions 
marked "overestimated" are regions where we have too much suppression due to the new form factors. 
These differences can be removed by modifying the non-Sudakov and Sudakov form factors. 



in Fig. [271 Of course by removing that contribution we get rid of the potentially large 
overestimate of the phase region. The missing parts would have had zero area if there 
were no recoils, since then Xi-i and Xi would have been equal. Thus we can fill in this 
region by associating a non-Sudakov ls!^J^ with the soft emissions as well. It is easy to 
see from Fig. [27] that the missed part extends down from -fCj-i to the angle of emission 
i — 1, i.e. to ^i-i- Thus for each soft emission we associate 

A;5"(i)=expf-5,r3£ r--'^) (B,4) 

which we have written in terms of the rescaled momenta as before. Again if Ki^i < zpi-i, 
(1B.4P will give an enhancement in some region but this is correct since in that case it will 
cancel an oversuppression coming from A^. 

From Fig. [27] we can see that there are still some regions in which we have an oversup- 
pression coming from both and from Ans associated with the hard emissions. Actually 
depending on the kinematics (when Ki < ZiPi), Ans will cancel the oversupression from 
As in some small region. This has been taken into account in Fig. [27] As one can see one 
overestimates the phase region when qi > Ki^i. Therefore in this case we can modify the 
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form factors as follows 



A^-'i(i) ^ expl^-asj^^^ \iq^>K,^i 

= exp -as \ — / (B.6) 

so that in each factor we integrate in energy up to yi instead of integrating up to Xi^\. In 
this way we minimize the overestimation of the true phase region. 

To summarize, the standard splitting probability for the emission of the ith gluon 
which is given by 



dP{i) = as dz,^ + ) As{mP^ - Zi^m-i) (B.7) 

irpf \ Zi I- Zi ' 



with As{i) given in (j2.23p and A„s(^) in (|2.2ip . is now replaced by 



dP{i) = as dzi — J 0{Ki^i - (1 - Zi)pi) \ + 

-n-pf \ \ Zi 1- Zi ' 



\ Zi i Zi 



where Af Anl\i), Aj«^'^(i) and A^f'^(i) are given by ([Ell), dEl]), (|R5]) and ([BiB]) 
respectively. 

It is really non-trivial to implement these changes in the type of numerical procedure 
that we use in this paper. However, they are suitable for implementing in a Monte Carlo 
procedure and as such it would be straightforward to implement them in CASCADE, and 
it would be interesting to see how large the effects are. One can also put in theta functions 
to make sure that for example z < 0.5 for hard emissions, and z > 0.5 for soft emissions. 

C. Re writ ing ^fig in CCFIVE ciiid the compfirison to BFKL 

In the literature one usually finds the statement that Eq. (|3.2p reduces to the BFKL 
equation in the high energy limit. In this appendix we would like to elaborate more on 
this. The key point to this statement is the derivation of A^s appropriate in the formal 
high energy limit, where energy is infinite, and importantly — > for fixed qq. 

Previously we derived A„<j following figure [25l In this case, however, we must consider 
what happens in case z is so small that the diagonal line through q, on which ^ is constant 
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Figure 28: The phase space for Ans in CCFM in case we take very large steps in z, such that z < qo/p. 



and equal to the angle of the real gluon q, cuts the line representing the soft cutoff qo 
before cutting the vertical line through k. This situation is illustrated in Fig. [28j The 
horizontal line at the bottom just indicates the momentum cut go- What happens here is 
that z < qo/p where p is the rescaled momentum of the real emitted gluon q. Actually 
since we now consider the small z limit, there is no need to make a distinction between 
the rescaled and the regular momenta, thus one might as well set p = q. From figure (o), 
one can see that A^^ is integrated over a region of total size 



A 



1 1, p\ 1 
In in — + In — :r in — 



qo 



In - In — — In^ — 



z 



ql 



(C.l) 



where in the last equality we have just set p = q. A„s is then given by exp(— agA). In the 
case shown in figure {h) we instead have 



A„, = exp(-as(-A + B + C)), 



(C.2) 



where 



A + B + C = Iniln^ + ln^ln— -In^ln^ 
z qo qo pqo qo Qq 

= In — In — In^ — 

z qo qo 



(C.3) 



where we have again set p = q in the last equality. We now remember that the virtual 
form factor is in BFKL given by ()2.6p . which equals 



/ - 1 1 , 

exp —Os in - in 

z qo 



(C.4) 
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Assuming now that the typical z are so small that the result in (|C.ip is valid, we can write 
the CCFM non-Sudakov form factor as 



A„,(z,A:,g) = A(z,A;) • A^(g) (C.5) 
where A(2;, A;) is the BFKL form factor, and 

A^(g) = expfa,ln2^) (C.6) 



90, 

with the " i?" standing for "Residual" . In this case Eq. (|3.2|) can be written as 

A{Y, k, q) =as r dy [ ^e{Y-y + Hq/q))A{Y - y, k)A^{q) A{y, k' , q), (C.7) 

Jo J T^q^ 

where A is the BFKL non-Sudakov form factor, while A-^(g) is given in ()C.6p and is 
independent of Y . Now of course if one considers the leading contributions in y — y, then 
the angular ordering theta function can be neglected and, since we moreover have — > 
one can let A^ ^ 1. In that case obviously the BFKL equation is recovered. 

It is obvious from the derivation here that the form for A„s is appropriate only for 
very high energies, and for ^ with fixed cutoff qq. Thus if one makes a certain 
calculation, for example jet production rates using q^ as a resolution scale, one needs to 
be careful considering the limit when gg ~^ 0. Obviously in that case we gradually go from 
the situation in Fig. [28] to the situation in Fig. [25j This means that we go from the result 
(ICT]) to the result in ([MI]), in which case the similarity to BFKL is somewhat lost. For 



similar discussions see also [38l . |47| . |39| . Of course in a sense this is a bit artificial since 

eik 



from the beginning we had that 5^g • S'j-, = A^^^^, but remember also the discussion at 



the end of section 



D. Increasing the accuracy of the real- virtual cancellations 

In this final appendix, we would like discuss a bit more carefully the definition of the 
/c_L -conserving emissions, and the real- virtual cancellations discussed in section [3l We were 
not very careful when defining the -conserving emissions, and especially the situation 
when the momenta of these emissions, q, become comparable to k was not discussed. To 
be more careful, one would have liked to define A;_|_-conserving emissions as emissions with 
say q < ak for some a which is strictly smaller than 1, as this would give a better accuracy. 
As the procedure in section [3] gives rise to an intercept higher than BFKL, one might try to 
thus increase the accuracy in the definition of the A;_|_-conserving emissions to derive a new 
equation. Indeed we saw in section [32] that a better agreement with BFKL was obtained 
in the asymptotic limit if a constant negative term proportional to the gluon distribution 
is added to the equation. As such a contribution is related to the virtual corrections, we 
might indeed guess that it has something to do with the real-virtual cancellations. 

The issue of improving the accuracy in the definition of the /cj^ -conserving emissions 
was discussed first in [s^, again considering only the possibility that k > k' . In our case 
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(a) (b) 

Figure 29: (a): Region D is the region where we have real k±-conserving emissions when the accuracy 
of the dehnition of k±-conserving emissions are improved. The shaded region is than the region left over 
by the real-virtual cancellation. Here q > k. It is seen that z' — k'^ /q^ and z" — ak/q. We have assumed 
z > z" , and the case z' < z" is discussed in the text, (b): Same as in figure (a) but now for q < k. 



we also look at the opposite case, again with the kinematical constraint included. The 
procedure is the same as before, and one can easily derive new evolution equations. As 
we later discuss, however, this procedure does not simplify the problem of solving the 
more exclusive equation in ()3.2p , although it sheds some light on the physical origin of the 
discrepancy between BFKL and the asymptotic behaviour of equations (j3.1ip and (j3.12|) . 

Since now the /c^-conserving emissions have momenta a times smaller than the cor- 
responding virtual propagator, they will not completely cancel An.s, but there will be a 
region left over from the cancellation. Assume first that k < q. Remember from Fig. [9] 
that Ans = exp{—asA) where A was defined in the figure. The region A from Fig. [9] 
is then the sum of the regions C and D in Fig. [29lfa). A = C + D. Region D is the 
region where we have the real A;_|_ -conserving emissions, and summing over these emissions 
we get a factor exp(Q;sZ)) and what is left from the multiplication with A„<j is therefore 
exp(— OsC). From Fig. [29lfa) we see that 
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lnlln-^-ln2A=lnilnl-ln2 4-. (D.l) 
zq^ aq z ak 



If instead k > q > ak, we see from Fig. [29{b) that region C is again given by the formula 
above. Thus this is the region left over by the cancellations. Actually, in Fig. [29lfa) we 
assumed that z' > z" which means that a < k/q. If instead a > k/q we get a slightly 
different situation, since the region C would then be replaced by a rectangular region. In 
that case it is easy to see that we get 



C" = lnlln-^.=lniln4T-ln2 4-+ln2-^ (D.2) 

zq'^ z a^ ak aq 
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Figure 30: The Mellin space eigenfunction of Eq. (|D.4|l for a — 0.5 (long dashed lines) and a = 0.1 (short 
dashed lines) for — > together with the BFKL eigenfunction (solid lines). 



which contains an extra term compared to the formula for C. Since in this case k < q, we 
can neglect the last term, however. Actually we could well approximate the regions C and 
C by the first term only (which contains ln(l/z)). This is indeed appropriate as we will 
below only study the asymptotic limit. We see that this term will generate an additional 
contribution in the right-hand side of the differential equation which reads 



Os In —J A{Y, k). 



(D.3) 



This is infact the term we were looking for in order to reduce the Mellin space eigenfunc- 
tion. However, note also that we have real emissions left over in the same region and these 
will modify the real kernel as well. We will write down the asymptotic equation below 
and study its Mellin transform. 

Before doing that, let us comment on the procedure of generating new equations. As 
an explicit constraint we can obviously use q > mm(^ak,ak'). The problem now, however, 
is that we can no longer justify the step from (j3.7p to (j3.8p where we dropped all the 
dependence on the third parameter in the argument of A. The dependence on q drops 
again on account of the kinematical constraint, but we can no longer guarantee that 
q > k' . If we nevertheless assume this then we lose the accuracy gained by introducing a. 
Infact if it was possible to drop all dependence on the third parameter by just using the 
kinematical constraint, then we would not need to bother cancelling A^^ since it would 
have been equally easy to solve (13. 2p directly. As far as the asymptotic behaviour is 
concerned, however, we are justified to go from (|3.7p to (|3.8p . (Unless we are looking at 
the — > limit, we would not be able to drop the second and third terms in (jP.ip and 
()D.2p anyway.) Thus to study the intercept we make the same steps as before and arrive 
at the asymptotic equation 



dvAiY, k) = as 



dk 



12 



1 



|A;2 - A:'2 



■ ha{K)A{Y, k') - a. In ^ A{Y, k), 



(D.4) 
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where now 

/,„(.) = 1 - - arctan (^^^^^_^-A^ j 0(2,/^ - (1 - - 1). (D.5) 

and K is defined as before. 

It is straightforward to study the Melhn space eigenfunction numerically. In Fig. [30] 
we plot the eigenfunction for a = 0.5 and a = 0.1 together with the BFKL result. As 
a decreases we indeed see that the curves move toward the BFKL eigenfunction. We 
basically get a constant shift around 7 = 0.5, and the behaviour at the endpoints 7 = 0, 1 
is not changed. For the intercept, one goes from 3.23 at a = 1 to 3.02 at a = 0.5, and 2.83 
at a = 0.1. The saddle point occurs at 7 = 0.5 in all cases. The saturation saddle point 
moves from 7^ 0.35 at a = 1 towards 7^ ~ 0.37 at a = 0.1. 
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